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FAILURE  RATE  ESTIMATION  USING  RANDOM  SMOOTHING 

by 

Douglas  R.  Miller 
Nozer  D.  Singpurwalla 


In  "this  note,  we ’investigate;  some  aspects  of  the  distribution 
theory  of  an  estimator  of  the  failure-rate  function.  We  clarify,  and  give 
a slightly  different  approach  to  some  asymptotic  results  of  Singpurwalla 
<1975)  for  the  failure  rate  process  from  exponential  lifetimes.  In  addi- 
tion, we  extend  these  results  to  the  case  of  general  lifetimes.  This 
leads  us  to  confidence  region^  for  "randomly  smoothed"  versions  of  the 
failure-rate  function,  for  both  finite  and  large  samples. 
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1.  Introduction 


Let  F be  the  distribution  function  of  a random  lifetime,  and 

let  r = - ~ log(l-F)  be  its  failure  rate  function.  The  purpose  of  this 

note  is  to  consider  estimators  of  r based  on  complete  and  censored  data, 
and  to  investigate  the  related  distribution  theory  for  both  the  finite  and 
the*  infinite  sample  cases. 

Consider  a collection  of  n identical  items  which  function  inde- 
pendently of  one  another.  Let  (L^)  be  the  time  to.  failure  (withdraw- 

al) of  the  ith  item.  We  assume  that  the  X^'s  are  independent  and  identi- 
cally distributed  with  distribution  function  F . The  L.'s  are  random 

withdrawal  times  which  may  have  any  joint  distribution,  but  must  be  inde- 
pendent of  the  X.’s  . Let  7 1 = minCX^,!.^)  ; thus,  if  7.  = X ^ , then 

the  ith  item  has  failed  at  X.  , and  if  Z = L.  , then  the  ith  item  was 

withdrawn  from  observation  at  L^  . Suppose  that  k failures  have  been 

observed  in  all,  and  let 


° = Z(0)  - Z(l)  - Z(2)  - Z(k) 


T-34  7 


be  the  ordered  failure  times.  Let  N (t)  be  the  total  number  of  items 

n 

t 

on  test  at  time  t , and  let  T (t)  = / N (u)du  be  the  "total  time  on 

0 n 

test"  at  time  t . 


We  shall  first  define  a "naive"  estimator  of  r , say  R , as 

n 


R"w  ' VZ(1>)  - Vz(i-1)>  • ziZw  (1.1) 


, z > Z 


00  ‘ 


Some  Comments  on  the  Naive  Estimator 

We  note  that  the  estimator  R is  the  reciprocal  of  the  total  time 

n 

on  test,  on  an  interval  over  which  it  is  defined.  The  estimator,  in  effect, 

assumes  a piece-wise  exponential  distribution  between  successive  failure 

points.  In  order  to  compute  R (z)  , we  must  look  "ahead"  to  the  failure 

n 

which  immediately  follows  z . Thus,  the  estimator  R (z)  is  a left- 

n 

continuous  step  function. 

Assuming  that  there  is  no  censoring  or  withdrawals  (i.e.,  when  k=n) 
and  assuming  that  n is  large,  Sethuraman  and  Singpurwalla  (1977)  show 
that  R (z)  is  not  consistent.  They  induce  consistency  by  smoothing  this 

il 

« * 

estimatoi  using  windows  of  fixed  width.  The  width  of  their  windows  depends, 
among  other  things,  on  the  sample  size  n . In  this  paper,  we  shall  consider 

another  smoothing  procedure.  However,  we  first  obtain  a right-continuous, 

! 

closed  form  version  of  our  estimator  (1.1). 
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be  the  r ight-continuous  inverse  of  F , with  F (0)  = 0 . Then,  for 
Z(i-1)  - Z < Z(i)  ' 

FnlFn(z)  = Z(i)  3nd  Fn1(Fn(z)  ' k_1)  = Z ( i- 1 ) * 

Thus,  given  k failures,  a right- continuous  version  of  our  naive  estimator, 
in  closer1  form,  is 


R (z)  = = , : 

T (F  F (z))  - T (F  (F  ( z ) - k i)) 
n n n n n n 


(1.2) 


An  analysis  of  the  right-  and  the  left-continuous  versions  of  R (z) 

n 

will  be  the  same,  and  thus  we  will  not  make  any  distinction  of  this  through- 
out our  paper.  However,  the  notation  in  Sections  3 and  4 will  be  simplified 
by  using  right-continuity. 

1.1  A Smoothed  Version  of  the  Naive  Estimator 

We  next  introduce  what  we  call  a "randomly  smoothed"  version  of  our 
naive  estimator.  Our  motivation  for  choosing  this  smoothing  technique  is 
given  in  Section  1.2. 

For  0 < h < 1 , let 


[hk(n)]  = mim{i,  'ah  integer  : i _>  hk(n)}  ; 

that  is,  [hk (n)  ] is  the  smallest  integer  greater  than  or  equal  to  hk(n) 

where  k(n)  , the  number  of  observed  failures,  is  random.  For  Z,  . < 

(l-l)  - 

z < Z . , we  define  a smoothed  estimator  R . (z)  as 
fi;  n,h 

a (7\  = fhk(n)] 

n,li  ' T (Z,..)  - T (Z,.  ....  • 

n (l)  n (i- [hk (n)  ] ) 

In  closed  form,  the  above  estimator  can  also  be  written  as 


R (2)  = [hk(n)  I 

n’h  Tn(FnlFn(z))  ~ C2)"'1)) 

n n n n n n 


(1.3) 


- 3 - 
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If  wo  set  j = [hk(n)]  , then  for  Z^.  ^ < z < Z^^  , 


R . (z)  = — r — _] Tj 

T (F  JF  (z))  - T (F  X(F  (z)  - jk"'  )) 


n n n 


n n n 


(3.4) 


In  the  above  expression,  j is  the  number  of  failure  intervals  used  in 
estimating  the  failure  rate  at  each  fixed  time  point. 

Wo  now  give  two  equivalent  expressions  for  R . in  terms  of  R , 

n , J n 

our  naive  estimator.  These  two  expressions  reveal  the  "smoothing"  more 
transparently.  In  Section  1.3  we  shall  show  that  for  Z^_^  £ z < Z^^  , 


and  that 


(R  .(z)f1  = \ Y -R-~ y 

n,j  1 R (Z..  \> 

J m=0  n (l-m) 


(1.5) 


R .(z) 
n,  J 


/ (l)  R (u)dT  (u) 
7 n n 

= (i-1) 

/ (l)  dT  (u) 

7 n 

(i-j> 


(1.6) 


Since  T (t)  is  the  total  time  on  test  at  time  t , Equation  (1.6) 

can  be  regarded  as  the  filtering  of  R (u)  through  a random  window  T (u)  , 

n n 

the  total  time  on  test  at  tim(e  u ..  We  shall  refer  to  this  as  the  "total 
time  on  test  random  window."  . In  the  following  section  we  present  consid- 
erations which  led  us  to  the  choice  of  our  smoothing  technique.  In  the 
sequel,  we  also  introduce  some  notation  which  will  be  used  later. 


1. 2 Justification  for  Smoothing  Technique 


In  retrospect,  we  state  that  a reasonable  justification  for  the 
choice  of  our  smoothing  technique  is  that  it  leads  us  to  a tractable 
distribution- free  theory.  However,  a primary  eonsiderat ion  which  led  us 
to  our  smoothing  technique  was  our  desire  to  clarify  and  to  expand  the 
asymptotic  results  of  Singpurwalla  (1975).  Of  particular  interest  is  his 


4 
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Theorem  4.2,  which  states  that  for  an  exponential  distribution,  the 
first  difference  of  the  reciprocal  of  the  estimated  failure  rate  converges 
weakly  to  a moving  average  process.  He  uses  this  theorem  to  give  a par- 
tial justification  for  an  empirical  (Box-Jenkins  type)  time  series  analysis 
of  the  estimated  failure  rate.  Singpurwalla  attempts  to  prove  this  theorem 
by  using  the  following  theorem  (see  Section  2): 

If  F is  exponential,  then  for  t E • • • »k(n)  , 

the  "total  time  on  test  process"  VM  . . . = T (Z..  . )/T  (Z,  . .)  converges 

[k(n)t ] n (l-l)  n k(n) 

weakly  to  the  Brownian  bridge  W°  [Billingsley  (1968)].  If  for  any  fixed 
h , 0 < h < 1 , we  define  the  difference  operator  * [0,1]  [ h , 1 ] ; 

V.  f (x)  = f(x)  - f(x-h)  , then  vf  W°  is  a moving  average  process.  Thus, 
h h 

considering  the  differences  of  the  total  time  on  test  process  V.,  , , , be- 

[k(n)t] 

comes  quite  natural. 

We  shall  see  in  Section  3 that  V,  V.,  . . , , h = j/n  , is  the  recip- 

h [k(n)t] 

rocal  of  our  randomly  smoothed  estimator  R .(z)  given  by  Equation  (1.4); 

n » J 

this  was  the  rationale  behind  the  choice  of  our  estimator. 

We  can  provide  some  further  insight  into  the  choice  of  our  smoothed 

estimator  if  we  note  that  the  right-hand  side  of  Equation  (1.4)  denotes 

the  total  number  of  failures  j > 1 divided  by  the  total  time  on  test 

( ’ 

between  Z...  and  Z,.  ..  . Thus,  our  smoothed  estimator  is  merely  an 

(i)  (i-j) 

extension  of  our  naive  estimator,  and  is  obtained  by  considering  a random 
interval  Z^  “ Z^  ^ , where  j = [hk(n)]  . 


1.3  Equivalence  Between  the  Versions  of  R . (z) 

n O 


It  now  remains  to  be  shown  that  (R  .(z))  equals  the  right-hand 

n * J 

side  of  Equation  (1.5),  and  that  Equations  (1.5)  and  (1.6)  are  indeed 
equival ent . 
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To  sec  that  the  former  is  true,  we  note  [cf.  Barlow  and  Campo 
(1975)],  that 


k( 


~r  (T  (Z,..)  - T (7.  ..))=/  (l)  F (u)du 

(n)  n (l)  n (l-j)  „ n 


'(i-j) 


and  that  R (z)  is  specified  by  Equation  (1.1). 
n 


In  order  to  verify  the  equivalence  between  Equations  (1.5)  and  (1.6), 
we  recall  a property  of  the  total  time  on  test.  This  property  [cf.  Barlow 
and  Campo  (1975)],  states  that  the  derivative  of  the  total  time  on  test  is 
the  reciprocal  of  the  failure  rate. 


1.4  Summary 


The  remainder  of  this  paper  is  devoted  to  investigating  the  distri- 
bution theory  of  R , and  demonstrating  its  usefulness, 
n , h 


In  Section 
we  apply  these  to 

expand  upon  some 
the  distribution 

present  a theory 
filtered  through 


2 we  present  some  preliminary  results,  and  in  Section  3 

R . when  F is  exponential.  Here  we  also  clarify  and 
n , n 

of  the  results  of  Singpurwalla.  In  Section  4,  we  consider 
theory  of  R^  ^ for  both  finite  and  large  samples,  and 

for  construction  confidence  bands  for  r , when  r is 
the  total  time  on  test  window. 

I 


In  Section  5,  we  discuss  the  calculation  of  critical  values  for 
use  in  constructing  confidence  bands.  In  Section  6 we  illustrate  our  tech- 
nique for  estimating  the  failure  rate  of  AC  generators,  based  on  some  real 
'life  data.  In  Section  7 we  generate  failure  data  from  a known  distribution 
and  use  it  to  estimate  the  failure  rate.  We  conclude  the  paper  with  some 
discussion  of  the  technique. 


2 . Preliminary  Distribution  Theory 

We  shall  assume  the  conditions  of  failure  and  random  withdrawals, 
as  stated  in  Section  1.  We  shall  also  adhere  to  the  notation  of  Section  1. 


6 
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The  text  of  our  paper  is  based  on  the  following  distribution-free 
result.  Lemma  2.1  is  a mild  extension  of  that  due  to  Barlow  and  Proschan 
(1969). 

Lemma  2.1:  For  any  distribution  F (F(0  )=0)  with  failure  rate 

r(*)  , the  random  variables 

Z . 

Y.  = / 1 r(u)N  (u)du  , i=l , 2, . . . , k(n) 

l 7 n 

( i- 1 ) 

— y 

are  independently  distributed  with  density  e . 


Proof:  The  proof  follows  that  of  Barlow  and  Proschan,  by  condition- 


ing over  the  values  of  L^  and  by  assuming  independence  between  the  L^ 
and  the  X.  . 

l 


// 


We  recall  that  k(n)  equals  the  observed  number  of  failures  from  a 
collection  of  n identical  items. 


Theorem  2.2:  When  n -*■  m , 

1 — - — If  (^(n)]  r(u)N  (u)du  - tk(n)l  , 0 < t < 1 1 {W(t)  , 0 < t < 1}  , 

| A(n)  \0  n 

where  ->  indicates  convergence  of  random  functions  and  W is  the  standard 
Brownian  motion  {Billingsley  (1968).]. 


■) 


Proof:  From  Lemma  (2.1), 

Z 


1 1 k (n) 


(tk(n) ] 

r(u)N  (ujdu  = £ Y.  ; 

n i=l  1 


(2.1) 


the  result  follows  from  Donsker’s  Theorem  (Billingsley,  p.  68). 


// 


Corollary  2.3:  When  F is  exponential  with  r(u)  = X , u > 0 , 


tl  , 0 < t < 1 


•+  { W(t ) , 0 < t < 1}  . 


7 
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For  completeness,  we  also  present  the  following  well-known  theorem 
[Pyke  (1969),  Shoraek  (1972)]. 


Theorem  2.4:  When  F is  exponential  with  r(u)  = X , u ^ 0 , 

£) 


L-T-y  /Tn(Z[tk(n)l)  \ 


< t < 1 } -►  {W  (t)  , 0 < t < 1}  , (2.2) 


where  W is  the  Brownian  bridge  [Billingsley  (1968)]. 


Proof : The  proof  follows  from  the  fact  that 


and  that 


[ tk (n)  ] 

W<„)]>  ‘ A b 


j k 

l Y / l Y , j=l , 2 k-1 

• _ 1 1 • _ 1 1 


i=l 


i=l 


are  distributed  as  the  (k-1)  order  statistics  from  a uniform  [0,1] 
distribution  [Karlin  (1966)].  These  in  turn  converge  to  a Brownian  bridge.// 

Theorem  2.5: 

Z 


(a) 
and 

(b) 


[ tn (k) ] 

/ (R  (u)  - r(u))N  (u)du  = 7 (1  - Y.) 

0 n n i=l  1 


[ tn(k) ] 


— / ftn(k) 1 (r  (u)  - r(u))N  ( u ) du  , 0 < t < 1 } + (W(t)  , 0 < t < l} 
A(n)  0 

Proof:  We  first  not'e  that 


1 


Z. .. 

/ 1 R (u) N ( u ) d u = -=—p 

hi-i)  " " (1>'  ‘"'“(‘-n'  Vn 


Z. 

x (y  \ T /7  x"  f 1 N (u)du  = 1 . 

T.O'/ix)  " T_(Z,,  ,N)  „ n 


Part  (a)  of  the  theorem  follows  from  this  and  Equation  (2.1).  Part  (b)  of 
the  theorem  follows  from  Donsker's  Theorem.  // 


8 


T-347 


In  terms  of  the  "total  time  on  test  process,"  part  (b)  of  the  above 
theorem  gives  us 

j—1—  f/[tk(n)1  R (u)dT  (u)  - r (u)dT  (u)l  , 0 < t < ij 

1 /k<iO  lo  n 0 " / “ " J 


(2.3) 


1 tk  (n)  ) ) JD 

-ZZZZ.  I (1  - Y.)  , 0 < t < 1}  ->  {W(t)  , 0 < t < 1} 

A(n)  i i-1  1 i ) 


3 . Asymptotic.  Theory  for  Exponential  Lifetimes 

In  this  section,  and  the  following  one,  we  derive  some  distribution 

theory  for  our  smoothed  estimator  R We  do  this  by  taking  differences  of 

n , h 

the  expressions  derived  in  Section  2. 

The  difference  operator  V,  was  defined  in  Section  1.2.  For  any 

h 

fixed  h , 0 < h < 1 , V is  a continuous  function  in  Skorohod's  D-topologv, 

h 

and  thus  when  V,  is  applied  to  the  statement  of  Corollary  2.3,  the  continu- 
h 

ous  mapping  theorem  (Billingsley,  Theorem  5.1)  applies  to  give 


(V,  AT  (7,  , , . .) 

“ V(nT~~~ 


- V 


) , h < t < 1 J ^ 


(V  W(t)  , h < t < 1} 


Since  V^t  = h , and  since  Equation  (1.3)  implies 

( 

n / 7 \ _ [hk(n)  ] 

Rh,n('(ttk(n)]-l))  ‘ V„<Zttk(n)])  " 

we  have  proved 


Theorem  3.1:  When  F is  exponential  with  r(u)  = A , u > 0 , 


A(n)  (rr1 — 

n,h 


(ftk(n)]-l) 


y - { , h < t < lj  - {Vh  -1-  W(t)  , h < t < 1}  , 


which  is  equivalent  to 


’ 

^ (rT(zA'7~7T  - i)  ’ 0 - ^ 1 " {vh  h w(t)  . £ t < 1}  . 

\ n,n  [ tk(n) ] / 
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l 


An  analogous  result  can  be  proved  by  taking  differences  of  (2.2). 
Singpurwalla  (1975)  does  this;  however,  he  uses  a sequence  of  difference 
operators  ‘>  0 • In  this  context.  Theorem  5.5  of  Billingsley  must 

be  used.  If  we  apply  the  operator  ^i/n  to  che  statement  of  Corollary 

2.3,  and  assume  that  k=n  , we  obtain 


[tn]} 


->  0 . 


Since  V,  , t 
l/n 


we  have 


l/n  , and  since  Equation  (1.1)  implies 


R (Z 


,)  = 


n^([tn]-l)'  Vl/nTn(Z[tn))  * 


1 

r R (zr„  ,) 

( n [tn] 


- < t < 1 
n — — 


JO 

->  0 . 


The  above  result  is  result  is  consistent  with  the  fact  that  the  naive 

estimator  R ^ has  a non-degenerate  limiting  distribution  fcf.  Sethuraman 

and  Singpurwalla  (1977)].  Similarly,  the  limit  in  Theorem  4.2  of  Singpur- 
walla should  be  0. 


3.1  Time  Series  Analysis  of  Failure  Rates 

The  value  of  performing  an  empirical  time  series  analysis  on  the 
estimated  failure  rate  was  demonstrated  by  Singpurwall a. (1975) . A goal 
there  was  to  obtain  a limiting  process  of  the  form:  a trend  plus  a Box- 

Jenkins  type  ARIMA  process.  For  complete  samples,  it  follows  from  Theorem 
3.1  that  1 

(R„,l.<Z[tn)))'1  * + ^ ; 

that  is,  a constant  plus  the  increments  of  a Brownian  motion,  which  is 
similar  to  discrete  white  noise.  An  analogous  result  is  observed  by 
Sethuraman  and  Singpurwalla  for  the  naive  estimator  R^  . Thus  it  appears 

that  the  desired  goal  of  obtaining  a limiting  process  of  the  form:  a trend 


r|S 


m 


v 
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plus  a Box-Jenklns  ARIMA  process  is  not  possible,  unless  some  additional 
structure  is  assumed. 


However,  we  conjecture  that  the  above  goal  could  possibly  be  realized 
if  the  sample  size  is  finite,  and  if  the  withdrawals  are  dependent  on  life- 
times. We  are  studying  this  possibility  in  the  light  of  the  following 
ramifications : 

1.  If  the  lifetimes  are  exponential,  they  do  not  age;  this 
makes  it  difficult  to  imagine  how  withdrawals  depend  on 
the  lifetime. 

2.  Without  making  restrictions  on  the  life  distributions, 
it  is  impossible  to  determine  from  the  data  whether  or 
not  withdrawals  are  dependent  on  lifetimes  [Miller 
(1976)]. 


4.  Distribution  Theory  for  General  Lifetimes 


As  in  Section  3,  we  apply  the  difference  operator  V to  Equation 

h 


(2.3),  and  noting  that  h = j(k(n))  , we  obtain 


f (ttk(n)l  Rn(u)dTn(u)  - / ([tk(n)1  r (u)dT^(u) 
Z([tk(n)]-j)  Z([tk(n)]-j) 


ltk(n)] 

l (i-V 

i=[tk(n) ]-j+l 

(4.1) 


We  next  define  a "randomly  smoothed"  version  of  r. 


Z, 


/ (l)  r<u)dT  (u) 

Z n 

(i.-i) z < z < z 

T (Z,^)  - T (Z,,_„)  * Z(i-j)  - Z(i)  • 


rT  (z)  = 

n.j  ^'“(i)'  ^'“(i-j) 


(4.2) 


I 


Combining  Equation  (4.2)  with  Equation  (1.6)  we  obtain,  for  ^ £ 

Z < Z(1)  , 


(Tn<7(0>-Tn(Z(l-j)>)(Rn,j(*>-'T  .<*>)'  l , > 

" * m=i-j+l 


n * J 


where  Y , m=l ,2, . . . ,k(n)  are  independent  exponential  random  variables 
m 

with  mean  1. 


11 
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The  analogous  result  for  R is 

n ) 


(t  (F-1F  (z))  - T (F  L(F  (z)-h))]  (r  (z) 
\nnn  nnn  / \ n,h 


n,  [hk(n) ] 


(z) 


[k(n)F  (z)] 
n 


m=[k(n)(F  (z)-h)]+l 
n 


[k(n)t] 

(l-y  ) = l (1-Y  ) 

m ra=[k(n)(t-h)]+l 


t=Fn(z) 


[k(n)t] 

= V l (1-Y  ) 
h m=0 


t=F  (z) 
n 


By  Donsker's  Theorem  and  the  continuous  mapping  theorem  we  obtain 


— (t  (F_1(t))  - T (F_1(t-h)))  (r  (F'^t))  - 
c"  \ n n n n / \ n,h  n 


A.(n) 


n,  [hk (n) ] 


<<'»)  • 


0 

h < t < 1 > -*• 


V,  W(t),  h < t < 1 . 

h — — 


(4.  4) 


Equation  (4.3)  can  be  used  to  find  the  finite  sample  confidence 


bounds  for  the  "randomly  smoothed"  failure  rate  function  r . The 


asymptotic  (n-*”)  confidence  bounds  for  r 


n,J 

can  be  obtained  from 


n [hk (n)  ] 


Equation  (4.4). 


Let  C+  , be  the  critical  value  such  that 
J ,k,a 


Pf  sup  ] (Y  -1)  < C+ 

( j<i<k(n)  m=  i- j+1  m.  J,k’ 


= 1-ci  , 


and  define  C.  , in  a similar  manner  for  tlie  infimum.  Then,  a 100(l-cx)% 

j , K , CX 


upper  confidence  bound  for  r^  is  obtained  from  Equation  (4.3)  as 


R.  ^ + c+  L.  - (T  (F  <z>>  ~ T (F_1(F  (z)-jk  S)j  1 > r (z)  . 

J,n  j ,k,a  \nn  n nnn  / — T . 


n,  J 


for  all  z,  7,  z < 7,,.  . 

(j-1)  - (k) 
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The  above  simplifies  to 


E°r  Z(j-1)  - Z < Z(k) 


R.  (z)(l  + C /j)  , (4.5) 

J » k , oc 

A lower  confidence  bound  is  similarly  defined. 


An  approximate  100(l-a)%  confidence  bound  is  obtained  from  Equation 
(4.4)  as 

Vn(*l  ♦ ATn7  < „i0  (t„<p;1<Fii<z)))  - Tn<Fn1<Fn(z)-h)))"1 

> r (z)  , 

n, [hk (n) ] 

-4* 

for  all  z , Z 


....  , < z < Z ,,  . , where  C,  is  the  critical 

([hk(n)J-l)  - (k)  h,°°,a 


value  such  that 


sup  V W(t)  < C } = 1-a 

h<t<l  h ' h> 


A lower  confidence  bound  and  a two-sided  confidence  bound  can  be  similarly 
obtained. 


5 . Critical  Values  for  Confidence  Bounds  for 
Randomly  Smoothed  Failure  Rate 

To  the  best  of  our  knowledge,  there  do  not  presently  exist  analytical 

methods  for  calculating  C.  , or  C*  . We  therefore  used  a Monte 

j,k,a  h,°°,a 

Carlo  method  to  calculate  selected  values  needed  in  Sections  6 and  7;  these 
are:  i 


1,55, .10  ~-9981 

C1 , 55 , . 10  ” 5,279 

1,55, .05  = ~'9990 

Cl,55,.05  = 6,008 

’ 5 , 5 5 , . 1 0 = ~4-06 

C5,55, .10=  7,887 

'5, 55, . 05  = ~4'21 

°5 , 55, . 05  = 9,02  * 
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The  values  are  based  on  one  run  consisting  of  20,000  replicates.  On  a 
second  run,  also  consisting  of  20,000  replicates,  it  was  determined  that 
the  four  two-sided  confidence  bands  using  the  above  values  would  have 
confidence  values  approximately  equal  to  .8035,  .8988,  .8094,  and  .9052, 
respectively.  Thus,  the  one-sided  critical  values  can  be  used  for  approx- 
imate two-sided  confidence  bands. 

6 . Application  to  Failure  Data  on  AC  Generators 

We  illustrate  our  estimation  and  smoothing  technique  by  considering 
some  failure  data  on  AC  generators  reported  in  NAILSC  Report  ILS  04-21-72. 

This  data  was  previously  considered  by  Castellino  and  Singpurwalla  (1973), 
and  by  Singpurwalla  (1975).  This  data  is  presented  in  Table  1 and  consists 
of  failure  and  removal  (withdrawal)  times. 

In  Figure  1 we  present  the  estimated  failure  rate  using  Equation 
(1.1).  The  confidence  bounds  using  Equation  (4.5)  are  so  wide  that  they 
are  not  useful.  For  instance,  a 90%  lower  bound  is  (1  - . 998l)R(z)  = 

.0019R(z)  , whereas  a 90%  upper  bound  is  (1  + 5.279)R(z)  = 6.279R(z)  . 

In  Figure  2 we  show  the  estimated  failure  rate  using  Equation  (1.4), 

R,.  . We  also  show  the  90%  upper  and  lower  confidence  bounds  for  r^  (•)  , 

’ n,5 

the  randomly  smoothed  version  of  the  true^ failure  rate.  We  remark  here  that 

the  estimator  R,.  was  calculated  using  the  current  interval  plus  the  two 

leading  and  the  two  lagging  intervals.  Clearly,  the  smoothing  leads  us  to 
narrower  bands  than  those  obtained  by  considering  single  intervals.  The  main 
issue  here  is  that  the  confidence  bands  pertain  to  a version  of  r , rather 

i 

, than  to  r itself.  Finally,  the  periodic  behavior  of  the  failure  rate, 

observed  by  Singpurwalla  (1975),  is  also  exhibited  by  Rr  . 

3 J 

7 . Application  to  Monte  Carlo  Example  and  a Discussion 

In  order  to  assess  the  performance  of  the  above  estimation  procedure 
in  a known  situation,  we  consider  failure  data  generated  from  the  failure 
rate  function 
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r (t)  = (l  + | + sin  2|)  / 12  . (7.1) 

In  Table  2 we  present  the  failure  times,  and  in  Figure  3 we  siiow  a 
plot  of  r(t)  . In  Figure  4 we  present  a plot  of  the  estimated  failure 
rate  R „ , together  with  the  90%  upper  and  lower  confidence  limits  for 

r_  . Tn  Figures  3 and  6,  we  show  the  "randomly  smoothed"  failure  rate 

n,5 

functions  r and  r , respectively.  If  we  superimpose  Figures  4 

n, 1 n,5 

and  6,  we  see  that  r falls  well  within  the  confidence  band. 

n,5 


7.1  Concluding  Remarks 

We  would  like  to  give  some  interpretation  to  the  randomly  smoothed 
failure  rate  r , smoothed  by  the  total  time  on  test  window. 

Suppose  that  we  first  consider  R^  our  naive  estimator,  which  gives 

a constant  value  over  failure  intervals.  Thus,  it  is  logical  to  consider 

this  as  an  estimator  of  r averaged  over  the  interval  of  interest;  however , 

this  is  essentially  r . Unless  some  regularity  conditions  are  imposed 

n,l 

on  r , it  will  be  impossible  to  obtain  confidence  bounds  for  r over  the 
interval.  Thus,  in  order  to  obtain  confidence  bounds,  we  consider  r^, 

< ' n«X 

instead  of  r . However,  the.  confidence  bound  on  . r is  so  wide  that 

n,l 

we  smooth  over  more  than  one  interval.  Thus  in  essence,  it  is  a small  step 
to  go  from  smoothing  over  o ine  interval  to  smoothing  over  several  intervals, 
i Finally,  even  though  the  smoothing  window  is  random  (it  depends  on  the  data), 
it  can  be  completely  specified  by  the  failure  data  and  is  thus  completely  known 
for  the  purpose  of  interpreting  the  confidence  band. 
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TABLE  1 

FAILURE  DATA  FOR  AC  GENERATORS 


Failure  Times 


1.0 

1.3 

3.0 

3.0 

3.8 

5.2 

24.6 

29.5 

30.6 

34.5 

38.9 

40.9 

70.0 

71.6 

100.5 

141.1 

160.8 

164.0 

200.3 

206.2 

212.5 

229.0 

252.8 

252.8 

274.3 

275.3 

282.2 

403.8 

435.3 

442.0 

466.6 

474.6 

495.0 

500.0 

502.5 

509.0 

520.7 

523.0 

599.0 

666.6 

676.0 

677.0 

744.8 

827.1 

852.0 

861.0 

950.4 

1097.3 

Censoring 

Times  (Withdrawal  Times) 

0.5 

0.5 

1.0 

1.3 

1.5 

2.6 

6.1 

6.2 

6.4 

6.5 

12.0 

12.7 

17.0 

17.2 

19.2 

19.5 

20.1 

23.7 

26.4 

28.6 

34.5 

36.6 

38.9 

43.4 

50.6 

53.7 

58.4 

59.3 

62.6 

64.9 

67.3 

69.3 

72.6 

72.6 

74.7 

74.8 

79.6 

80.0 

81.0 

83.4 

83.4 

84.3 

89.0 

90.1 

90.5 

91.3 

92.8 

100.2 

107.1 

109.0 

111.8 

112.1 

113.8 

117.0 

126.6 

131.2 

131.5 

132.0 

134.3 

135.0 

142.1 

149.2 

149.7 

157.1 

158.4 

161.0 

164.0 

166.4 

173.5 

174.0 

184.4 

187.2 

211.3 

216.5 

218.2 

225.3 

228.9 

233.3 

243.1 

261.0 

265.3 

265.7 

268.0 

268.6 

272.7 

274.6 

274.9 

, 275.2 

283.5 

292.0 

306.2 

314.0 

319.0 

320.7 

321.4 

328.7 

338.2 

347.8 

349.0 

355.1 

359.8 

371.8 

381.0 

382.8 

389.0 

393.2 

402.6 

403.7 

406.8 

411.3 

414.2 

428.8 

435.3 

441.0 

461.1 

465.7 

466.6 

479 ; 6 

498.8 

505.8 

512.0 

514.8 

518. ‘3 

585.2 

605.2 

606.4 

627.7 

627.7 

630.0 

645.9 

674.4 

684.3 

689.9 

701.7 

704.4 

724.4 

732.7 

752.5 

784.8 

784.9 

824.4 

864.3 

881.7 

898.7 

957.0 

983.1 

1064.1 

1147.3 

1173.6  1762.2 
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10.6 

43.2 

167.0 

272.7 

444.0 

510.8 
703.6 


5.2 

16.2 

24.7 

45.6 

65.9 
79.2 

86.9 
105 . ^ 
125.9 

140.3 

164.0 
191.6 

237.4 

270.5 

301.8 

330.6 

373.6 

406.1 

446.1 

510.8 
617.4 

685.2 

771.2 
952.0 
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TABLE  2 

MONTE  CARLO  FAILURE  TIMES 


.127 

.406 

.425 

.635 

.669 

1.041 

1.132 

1.156 

1.488 

1.788 

1.907 

1.991 

2.056 

2.122 

2.362 

2.366 

2.473 

2.604 

2.900 

2.970 

3.083 

3.372 

3.557 

3.943 

4.012 

4.352 

5.103 

5.615 

6.843 

6.928 

7.280 

7.942 

8.310 

8.311 

8.546 

8.792 

8.930 

9.137 

9.241 

9.561 

9.904 

9.988 

10.375 

11.507 

11.735 

12.748 

13.301 

13.508 

14.435 

15.518 

16.857 

17.266 

17.426 

17.975 

18.259 

i 
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